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Scattering rates for a Goldreich-Sridhar (GS) spectrum of anisotropic, incompressible, magneto- 
hydrodynamic turbulence are calculated in the quasilinear approximation. Because the small-scale 
fluctuations are constrained to have wave vectors nearly perpendicular to the background mag- 
netic field, scattering is too weak to provide either the mean free paths commonly used in Galactic 
cosmic-ray propagation models or the mean free paths required for acceleration of cosmic rays at 
quasi-parallel shocks. Where strong pitch-angle scattering occurs, it is due to fluctuations not de- 
scribed by the GS spectrum, such as fluctuations generated by streaming cosmic rays. 



The scattering of energetic particles by turbulent magnetic and electric fields plays an important role in the accel- 
eration and propagation of cosmic rays Jll^ . The turbulent fields responsible for cosmic-ray scattering can be excited 
by the cosmic rays themselves or by some mechanism that is independent of the cosmic rays. This paper focuses 
upon the latter case. In previous treatments of scattering, different turbulence models have been used, including 
fluctuations with wave vectors k parallel to the ambient large-scale magnetic fleld Bq (slab symmetry) or perpen- 
dicular to Bq (2D), or power spectra that are isotropic in fc-space ifZHlOll. On the other hand, a number of studies 
suggest that in magnetohydrodynamic (MHD) turbulence excited by large-scale stirring, small-scale fluctuations have 
non-zero values of fc|| that are <^ fcj^, where fcx and fc|| are the components of k _L and || to Bq |12-1J]. In this paper, 
the quasilinear approximation |11| is used to calculate general scattering rates for incompressible MHD turbulence 
and also shear-Alfvenic turbulence on the non-MHD scales shorter than the collisional mean free path of thermal 
particles These rates are then evaluated for the Goldreich-Sridhar power spectrum |jl^, which has significant 

power at small scales only for k± ^ k\\. The condition k± ^ k\\ is found to significantly decrease the efficiency of 
pitch-angle scattering relative to the slab-symmetric and isotropic cases. Astrophysical applications and limitations 
of quasilinear theory (QLT) are discussed. 

It is assumed that there is an inertial-range spectrum of fluctuations extending from some large scale I to a much 
smaller scale d, with the fluctuations at scales I dominating the total magnetic energy. Only cosmic rays with 
gyroradii p ^ I are considered. A scale V is introduced, with p <^ I' <^l. The energetically dominant fluctuations on 
scales > /' are treated as a uniform field Bq. The magnetic fluctuations on scales < I' , denoted Bi, are small compared 
to Bq and are treated using QLT. It can be verified a posteriori that the QLT scattering rates are independent of I' 
to lowest order in p/l. In contrast to most previous treatments, the turbulence is treated as strong, in the sense that 
fluctuations decorrelate in one linear wave period. 

In QLT, the turbulence causes the cosmic rays to diffuse in momentum space, with the diffusion coefficients deter- 
mined by the statistical properties of the turbulence |15[, 
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where / is the cosmic-ray distribution function averaged over the small scales of the fluctuating fields, p is momentum, 
9 (the pitch angle) is the angle between p and Bq, and ^ = cos 6*. In equation (|I|) it has been assumed that the length 
scale characterizing variations in / is large compared to p, so that / can be taken to be independent of gyrophase. 
At each k, Bi and the incompressible turbulent velocity Ui are decomposed into shear- Alfven and pseudo-Alfven 
components by projecting along the appropriate polarization vectors ]T^ . These components are denoted respectively 
by the superscripts s and p, so that 

Bi(k,t) = B^i(k,t)+BP(k,t), (2) 

with an analogous equation for Ui(k, t). The electric field is given by Ohm's Law, Ei = — (I/c)Ui x Bq. The 
normalized power spectra of the shear- Alfven modes are given by 
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with analogous equations for the pseudo-Alfven modes, where va is the Alfven speed associated with Bo, and (. . .} 
denotes an ensemble average. It is assumed that the turbulence is homogeneous and stationary, that (Bi(x, i)Bi(x + 
r, t+r)) — (Bi(x, t)Bi(x— r, t+r)) with analogous equations for (UiUi) and (UiBi) (no magnetic or kinetic hehcity), 
that (Ui(x, t)Bi(x + r,t + r)} = (Ui(x, t)Bi(x + r,t — r)) (which gives D^p — Dp^), and that the shear-Alfven and 
pseudo-Alfven modes are statistically independent. 

The contributions to the momentum diffusion coefficients from the shear-Alfven modes and pseudo-Alfven modes 
are, respectively, 
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where S = va/v, z — k±p, p — vj_/Q, is the cosmic-ray gyrofrequency, vj_ and W|| are the cosmic-ray velocity 
components -L and || to Bq, L is the dimension of a window function that multiplies the variables before a Fourier 
transform is taken, and the arguments of K, C, and M are (k, r). Since the shear-Alfven and pseudo-Alfven modes 
are statistically independent, D^^ = D^^ + D^^, etc. Equations (^ and (0) are derived using a standard method based 
on the linearized Vlasov equation jll], modified to treat strongly turbulent fluctuations instead of waves satisfying 
linear dispersion relations. Alternatively, they can be derived from equations (7a), (7b), and (7c) of |15|, if one notes 
the typographical error on the eighth line of equation (7a), namely, that Qr\\ should instead be Q||/;. 
A Goldreich-Sridhar spectrum of strong, anisotropic MHD turbulence is now assumed, with 
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for {I') ^ < k± < d ^ with d — > 0, where Tk = {l/vA){k±l) "^^^ is the Lagrangian correlation time appropriate for 
strong anisotropic incompressible MHD turbulence, and 



if \x\ < 1 
if \x\ > 1 



(9) 



The spectrum of equation (jg) is also taken to describe the fluctuations on scales between V and I, and the normalization 
has been chosen so that the total magnetic energy J^^^ k±dk± dk^^M^{k±, k\i ,0)Bq/4: = L^Bq/Stt. At small scales 
in MHD turbulence, there is equipartition between magnetic and kinetic energies, so that = and = M^. 
It is assumed that AfP = ]\/P and 0^ = 0^ — aM^, where the arguments of each of these spectra are (A:_l, fc|| , t), and 
where the fractional cross helicity a G (—1,1) is independent of k. 
In many applications, there are two small parameters. 



For sin 6* > e^/^, 
order in e and (5, 
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one finds from equations m) and (m) and the assumed forms of the power spectra that to lowest 
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The terms on the right-hand sides of equations (^l]) and ( [T^ proportional to e^/^ correspond to fluctuations satisfying 
the magnetostatic gyroresonance condition fc||z;|| = nil, which states that the Doppler-shifted frequency of a static 
magnetic fluctuation in the reference frame of an energetic particle's motion along Bq is an integral multiple of the 
particle's gyrofrequency. A fluctuation is seen as static when a cosmic ray passes through one wavelength of the 
fluctuation in a time (fc||i;||)~^ ^ Tk. The gyroresonant terms in equations dll] ) and ( |l2| ) arc much smaller than in the 

case of slab-symmetric or isotropic turbulence for sin 6* ^ e^/^ because equations (^) and (H) imply that k± > /cjj^^Z^/^, 

so that fluctuations satisfying = nfl also satisfy k^p > n^^'^e^^^^ sin9 cos"^^^ 9 ^ 1. The condition kj^p ^ 1 
implies that a cosmic ray traverses many uncorrelated fluctuations of the required fc|| during a single gyro orbit. 
The effects of these uncorrelated fluctuations tend to cancel. The weakening of gyroresonant scattering due to this 
gyro-orbit averaging would occur for any power spectrum in which all fluctuations on scales <C I satisfy k± ^ . 
The terms on the right-hand sides of equations ( pi] ) and (|l^) proportional to (— Ine)^ correspond to non-resonant 
interactions. In equation ( [T^ ) the non-resonant term arises from the n = term in equation which represents the 
effects of the magnetic- mirror force of the pseudo- Alfven modes (transit-time damping) . This term becomes large as 
9 n/2, since as Dy va particles can "surf" magnetic mirrors moving at speeds ~ va more effectively. 

For sin 6* -C e^/^, scattering is dominated by magnetostatic gyroresonant interactions with shear- Alfven modes, and 
to lowest order in e and 6 

\ /.X n9^ 

pp 

Although D^^ vanishes as 6* — > 0, the pitch-angle scattering frequency i/ = 2D^^/{1 — ^^) (which unlike Z?^^ is 
independent of 6 for isotropic scattering) approaches {tt/ A){v/l) as — > 0. Gyroresonant interactions are stronger for 
9 < e^/^ than for 9 S> e^^^, because when 9 < e^^^ modes satisfying fc||ti|| = nfl also satisfy k±p ^ 1, so that a cosmic 
ray doesn't traverse many uncorrelated resonant modes during a single gyro orbit. 
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FIG. 1. The e dependence of the pitch-angle scattering frequency i/ = 2(D|j -I- -D|'^)/(l - C^) for e = 5 = 10"^. The xs 
indicate numerical evaluations of equations and (j^, the solid line gives the analytic results of equations (^) and (|l^, and 
the dashed lines give the limiting value from equation (^3|) of u as sin^ 0. 

In figure |l], the pitch-angle scattering frequency from equations ( |ll|) and ( |l^ ) is plotted with the solid line, and the 
limiting value of i' as sin 9 — > from equation (O) is given by the dashed line. The xs indicate numerical evaluations 
of ly from equations (^) and (^ for the assumed spectra in which only those terms in the infinite sum with \n\ < 10 
are kept and in which I' — Q.ll (see introduction). The values e = 10~^ and 6 = 10~^ have been used, (ct, which only 
weakly affects v, has been set to 0.) The characteristic values v ~ v/l for sin 6' < e^/^, ly ^ S~^v/l for Itt — 6*1 < S, and 
v ^ [(— In e)S + e^/'^]v/l for moderate pitch angles can be extrapolated to all values of e and 5 much less than 1 in the 
quasilinear approximation. 
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When V is sufBciently large, /n becomes nearly isotropic, and the pitch-angle-averaged distribution / can be treated 
in the diffusion approximation |7|,|l6|. 



dt ~ di r" 91 )^ p^dp "dp 



(14) 

where / is distance along a field line, and the ellipsis indicates the omission of the advection and adiabatic-acceleration 



terms. The coefficient of spatial diffusion along Bq is given by |16 
To lowest order in e and S when e^^^ <C (— lne)(5, QLT gives 
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[In this paper, from equation ([T^) on, MP = IVP; however, if either MP or is set to zero, equation (16) becomes 
n'^^'^ = vl{—2S\ii€)~^.] To lowest order in e and S when (— \ne)6 <^ e"!"^ QLT gives 

ci = (7r/22)csc(57!^l)65/"[(2/13)y;,';°^in-^/2]-6/ii ^ 9.88. The pitch-angle-averaged momentum diffusion coeffi- 
cient in equation (HJ) is given by p6| 
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To lowest order in e and (5, QLT gives 



^ = (V24)(l -a2)(-lne)p^«i/(«0. (19) 

When — 1, Dp vanishes since the small-scale fluctuations all travel in a single direction along Bq at the speed va, 
and, in the reference frame that follows their motion, particle energies are conserved. 

Although QLT is a useful and standard tool, it suffers from important inaccuracies. QLT assumes that during 
the time a particle is correlated with a turbulent fluctuation, the orbit of that particle is the same as in a uniform 
magnetic field. However, field-strength fluctuations A|i3| with A|i3|/|i?| = a <C 1 magnetically trap cosmic rays 
with 1^1 ^ a^l"^ . (For incompressible turbulent fluctuations, which have phase velocities ~ va along the magnetic 
fleld, and for cosmic rays with v vai the trapping condition is essentially the same as if the fluctuations were 
stationary.) The trajectories of such trapped particles differ greatly from the trajectories of particles in a uniform 
fleld, violating the QLT assumptions. Because the integral in equation ( |l^ ) is dominated by values of |^| ^ (5 ^ 1 
for which trapping is important, the value of in equation ( p^ is unreliable. Similarly, when (— lne)(5 <g; e^^^, 
the integral in equation (|l^) is dominated by small |^|, and thus the value of K|| in equation (ll^ is unreliable. In 
addition, assuming an unperturbed particle orbit in the presence of a slowly and non-periodically varying Ei or Bi 
leads to spurious changes in a particle's magnetic moment — mv\l2B^^ which, as an adiabatic invariant, should 
be virtually conserved when Ei and Bi vary on a time scale S> The non-resonant terms in equations ( pi] ) and 

( |l^ ) arise from slowly varying modes and imply such spurious changes in /i, thereby significantly overestimating non- 
resonant pitch-angle scattering. Since K|| in equation (|l^) is determined by this non-resonant pitch-angle scattering, 
equation ( p^ underestimates 

Although the QLT results for the key particle-transport coefficients are inaccurate, QLT does show that resonant 
scattering by MHD turbulence with fc_L 3> fc|| is much weaker than resonant scattering by slab-symmetric or isotropic 
fluctuations. Moreover, equation ( p^ ) as a lower bound on has important implications. If — 5/iG, I — 100 pc, and 
VA = 10^ cm/s [parameters characteristic of the interstellar medium (ISM)], then e = 2.2 x 10~^i?GeV for a relativistic 
proton, where E'goV is the proton's energy in GeV, and (5 = 3.3 x 10~^. If E'goV ^ 10^, then e^/^ -C (— Ine)^, and 
equation ( |T^ ) gives a lower limit to the scattering mean free path of 430 kpc x (20 — ln£^Gov)~^- This value 

is so large that if the power-law spectrum of interstellar turbulence inferred from observations is described by 
equation (@), then some mechanism besides such turbulence must be invoked to explain the confinement of cosmic 



4 



rays to the Galaxy At energies < 10^ — 10^ GeV, such a mechanism is provided by resonant MHD waves that 
cosmic rays themselves excite, but at higher energies it is believed that self-confinement does not work H,^. For 
-E-GcV > lO'^ — 10'^, cosmic-ray confinement and isotropization can be explained even if turbulent scattering is weak if 



one takes into account molecular-cloud magnetic mirrors [|18[. If the interstellar turbulence generated by supernovae 
and stellar winds is described by equation (ph, the inefficient scattering associated with such turbulence may indicate 
that quasi-parallel shocks are unable to accelerate cosmic rays up to the ~ 10^ GeV energies at the "knee" of the 
galactic cosmic-ray energy spectrum P,[l9t, although this suggestion is controversial |20|. Quasi-perpendicular shocks, 



however, may be able to accelerate cosmic rays to the knee and beyond pi 22 
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